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INTRODUCTION 


OVERVIEW  -  ACOUSTIC  INSTABILITY 

Oscillatory  flows  in  ducts  can  be  sustained  by  a  variety  of  interactions.  These 
range  from  purely  fluid-dynamically  induced  motions  (e.g.,  large  scale  vortex 
shedding  or  shockwave  instability),  to  the  more  energetic  motions  possible  with 
diabatic  flows.  Particularly  severe  oscillations  can  occur  in  solid  propellant  rocket 
chambers,  for  example,  wherein  the  oscillatory  motion  can  be  driven  by  interactions 
with  the  substantial  energy  release  inherent  in  near-surface  combustion  processes. 

Prior  analytical  work  in  combustion  instability  in  solid  rockets  has  identified 
some  of  the  mechanisms  which  can  produce  acoustic  velocity  coupling  to  the  overall 
instability  process  (see,  for  example,  Culick1-2).  These  studies  have  shown  that  for  a 
simple  longitudinal  standing  aeon  ic  wave  in  a  duct,  a  Rayleigh  diabatic  instability 
criterion  can  result  which  is  dependent  upon  the  oscillatory  motion  of  the  gas 
column.  Since  the  energy  release  occurs  in  the  near  surface  region,  analytical  work 
has  addressed  acoustic  boundary-layer  effects  on  propellant  combustion  response 
(see,  for  example,  Lengele3).  Implicit  in  several  response  function  analyses  is  the 
assumption  that  the  acoustic  boundary-layers  behave  quasi-steadily  and  in  phase 
with  the  longitudinal  acoustic  velocity  outside  the  boundary  layer. 

Recent  work  on  combustion-flowfield  interactions  in  solid  rocket  chambers 
has  analytically  and  experimentally  examined  fundamental  fluid-dynamic  aspects 
of  mean  flow,  acoustic  wave  and  turbulence  behavior.  Hydrodynamically  modeling 
the  flow  with  a  semi-enclosed,  porous-walled  duct  (with  large  injection  through  the 
wall),  Brown  et  al4  experimentally  confirmed  the  transitional  behavior  of  the  mean 
flow  predicted  by  Beddini5.  The  major  emphasis  of  this  non-reactive  flow 
experiment  was  to  investigate  the  effect  of  low  and  moderate  amplitude  forced 
acoustic  oscillations  on  tin  flow,  using  both  hot-wire  anemometry  and  surface- 
mounted  hot-film  sensors. 

Basing  their  conclusions  on  the  surface  sensor  measurements,  Brown,  et  al  4 
suggested  that  near  surface  turbulence  produced  by  the  mean-flow  transition 
process  appeared  to  “destroy”  the  coherent  response  of  the  sensors  downstream  of 
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nsition.  Consequently,  for  low  acoustic  amplitudes,  a  surface  response  capable  of 
luring  instability  would  most  probably  originate  in  the  head  end  (pre-transition) 
:tion  of  a  rocket  chamber.  However,  it  is  noted  that  for  large  axial  distances,  low 
iplitude  acoustic  signals  appear  to  be  resurrected  in  the  post  transition  region, 
wer  amplitude  harmonics  of  the  driving  frequency  were  observed  by  Brown  in 
di  the  surface  and  hot  wire  measurements.  Their  existence  is  consistent  with  the 
nlinear  behavior  of  a  (potentially  thick)  acoustic  boundary  layer  (Stokes  layer). 

Analyses  of  laminar  acoustic  boundary-layer  phenomena  involved  in 
apellant  response  have  been  offered  by  Flandro6;  Glick  and  Renie7;  Ben  Reuven8 
d  Hedge,  et  al.9.  Zinn  and  colleagues10  have  obtained  Schlieren  and  radiation 
lasurements  for  a  premixed  reactive  Stokes  Layer  on  a  porous  plate.  Baum  and 
wine11  utilized  a  full  unsteady  Navier-Stokes  solution  method;  their  results 
owed  a  Stokes  layer  persisting  even  in  the  presence  of  strong  injection  velocities, 
le  thickness  of  the  Stokes  layer  was  of  the  same  order  as  the  Stokes  layer  for 
ninjected  flow.  The  presence  of  vortices  above  the  acoustic  boundary  layer  were 
so  indicated  in  the  results.  Although  several  aspects  of  these  analyses  are 
teresting,  the  following  important  point  is  noted.  Simple  scaling  estimates  of  the 
minar  acoustic  boundary  layer  height  (the  Stokesian  thickness)  8a  =  "s/  v  /  f 
dicate  that  8a  is  well  above  the  gas-phase  flame  height  except  for  very  high 
?quencies  (f>  10  Hz)  .  Consequently,  minimal  acoustic  velocity  interaction  with 
mbustion  would  be  expected  at  low  to  intermediate  frequencies,  as  indicated,  for 
ample,  by  the  more  detailed  analysis  of  Flandro6.  Although  this  does  not  preclude 
locity-coupled  instability  resulting  from  laminar  interactions  (whether  linear  or 
inlinear),  it  does  suggest  that  other  mechanisms  of  interaction  should  be  explored. 

It  is  well  established  that  for  piston  driven  closed  duct  flow's  (viz.,  in  which 
ere  is  negligible  mean  axial  flow  *,  longitudinal  acoustic  waves  of  a  few  percent 
lative  pressure  amplitude  can  induce  turbulence  within  the  Stokes  layer.  The 
erature  review  and  independent  data  of  Merkli  and  Thomann12,  for  example, 
nfirm  a  critical  amplitude  which  varies  as  Vf  for  this  type  of  simple  acoustic 
otion.  Thus,  for  fixed  amplitude,  lower  frequency  disturbances  are  more  likely  to 
oduce  turbulence. 

It  is  known  from  studies  of  quasi-steady  solid  propellant/flowfield 
teractions  that  even  small  levels  of  turbulence  within  the  combustion  zone  can 
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appreciably  enhance  propellant  combustion  rates,  predominantly  by  increasing  heat 
transfer  from  the  flame  zone  back  to  the  propellant  surface.  Whether  acoustically 
induced  turbularization  can  occur  in  actual  rocket  chambers  depends  upon  several 
Complex  effects  including  the  presence  of  large  injection  rates,  surface  roughness, 
and  of  course,  the  specific  type  of  combustion  process.  In  one  experiment13,  values 
of  the  threshold  velocity  for  propellant  response  (acoustic  erosivity)  were  lowered  by 
a  factor  of  two  relative  to  steady  state  conditions.  Acoustic  turbularization  of  the 
propellant  combustion  zone  was  also  suspected  by  Medvedev  and  Revyagin14,  based 
on  their  T-burner  data  and  the  data  of  (Vice13. 

OVERVIEW  -  ACOUSTIC  EROSION 

Solid  rocket  motors  are  susceptible  to  several  types  of  flowfield  related 
combustion  anomalies.  SleaJv-state  erosion  and  acoustic  coupling  are  examples  of 
combustion  enhancements  wine  can  lead  to  anomalous  burning  or,  in  some 
instances,  catastrophic  failure  o!  the  motor  due  to  increased  chamber  pressures. 
The  steady-state  erosive  burning  problem  (also  referred  to  herein  as  erosive 
combustion  or  steady-state  erosion)  has  been  addressed  in  past  studies  and  is 
generally  believed  to  be  a  consequence  of  the  turbulent  mean  flowfield  enhancing 
the  combi  tion  an  1  heat  transfer  processes  by  increasing  the  effective  thermal 
diffusivity16.  Steady-state  erosive  burning  is  therefore  most  likely  to  occur  in  rocket 
motors  having  large  length  to  diameter  ratios  in  which  the  flowfields  attain  large 
velocities.  It  has  also  been  known  for  some  time  that  erosive  combustion  can  occur 
without  the  presence  of  an  appreciable  mean  flowfield,  specifically,  in  acoustic 
environments  found  in  the  center-vented  T-burner  experiments  of  Medvedev  and 
Revyagin14  and  Crump  and  Price1*1-17 

RESEARCH  OBJECTIVES 

The  objective  of  the  present  analysis  is  to  consider,  under  several  simplifying 
assumptions,  the  general  trends  of  acoustic  boundary-layer  turbularization  on  a 
permeable  surface  in  the  presence  of  large  injection  rates.  'The  fundamental 
nonreactive  problem  is  considered  first.  'The  analysis  is  then  extended  to  include 
the  unsteady  heat  conduction  and  combustion  of  the  condensed  phase.  The 
phenomena  of  interest  are  schematically  depicted  in  Figure  1  .  General  trends  of 


reactive  acoustic  boundary-layer  turbularization  on  a  propellant  surface  under 
nominal  homogenous  propellant  conditions  are  investigated  together  with  the 
resulting  response.  * 


Mean  Flow 
Profiles 


R  .  M  .  S  .  Acoustic 
Profiles 


Axial  \ 
Velocity  \ 

Turbulence 

Intensity 

V 

Tempei 

nature 

i 

i 

m 

t” 

< 

1 

■ 

<  T  > 

i 

i 

■ 

W//, 

Ts  ' 

rf 

wm 

T„ 


Figure  1.  Mean  flow  and  acoustic  profiles  at  a  pressure  node. 


The  combined  response  of  solid  propellant  combustion,  due  to  both  acoustic 
pressure  oscillations  and  the  acompanying  turbulent  flow,  is  also  investigated. 
Under  several  simplifying  assumptions,  the  general  trends  of  the  combined 
response  of  the  solid  propellant  at  various  positions  within  a  center  vented  T- 
burner,  or  low  Mach  number  region  of  a  solid  rocket,  are  determined.  A  comparison 
of  acoustic  and  steady-state  erosive  combustion  response  is  also  performed,  since  a 
close  relation  between  these  phenomena  has  been  assumed  in  prior  studies2. 
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ANALYSIS 


DECOMPOSITION  OF  DEPENDENT  VARIABLES 


The  motion  of  a  perfect  gas  described  by  the  Navier-Stokes  equations  is 
considered.  An  arbitrary  dependent  variable,  g(xpt),  is  decomposed  according  to  the 
notation 

g  (x,  ,  t)  =g(\,  ,  t)  +  g'(  x,  ,  t) 

with 

g  =  <g<  x,  ■  t)>  +  g"(  x,  ,  t). 


In  these  equations,  g  is  the  turbulent  fluctuation,  g  is  the  ensemble  average,  g"  is 
the  acoustic  (deterministic)  component  and  <g>  is  the  long-time  mean.  The 
remainder  of  this  section  is  diva;  d  into  two  parts,  the  first  considering  order-of- 
magnitude  analysis  of  the  problem  and  derivation  of  an  approximate  relation  for 
critical  amplitude  of  transition  from  the  turbulence  model  equations.  The  second 
part  discusses  the  equations  and  method  used  for  a  more  comprehensive  numerical 
solution  of  the  acoustic  boundary-layer  in  the  presence  of  transpiration. 

ORDER-OF-MAGNITUDE  ANALYSIS  18 

Neglecting  third  order  correlations  and  the  axial  derivatives  of  molecular  and 
turbulent  stresses,  the  <  >  average  of  the  axial  momentum  equation  is 


—  (<pxuxu>+<pxu"u">+<uxp"u '>+<uxp"u">) 
d\ 


+ —  (<pxuxv>+<pxu"v  ">+<uxp”v">+<vxp  "u">) 
dy 


d<p>  ^  d  (<g>  d  <u»  d<T> 
fix  fly  f)y  dy 


(1) 
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where  the  mean  turbulent  shear  stress  is  <t>  =  -<pxu'v'>  (neglecting  turbulent 
density  correlations).  The  equation  for  the  acoustic  component  (u  ")  is,  to  first  order: 


9u"  „  ,  „  d 

<P>— — +  (<P>U  +  <u>p  )+U  - (<pxu>) 

0t  0X 


0u"  „  „.0<u> 

+  <pxu> - +  (<p>v  +  <v>p  ) - 

dx  dy 


„d  ,  „  du"  dp"  d  0u"  dx" 

+  u  — (<pxv>)+  <pxv>- —  = - +  — (<p> — — )  +  —  (2) 

dy  dy  dx  dy  dy  dy 


where  t"  =  -<p>u'  v' . 

In  equation  (1),  the  second-order  correlations  arising  from  the  unsteady 
motion  can  produce  virtual  stresses  and  convective  fluxes  analogous  to  those 
produced  by  the  turbulent  motion  alone.  These  correlations  induce  the  classical 
phenomenon  of  “acoustic  streaming”,19-20  wherein  a  component  of  the  mean  flow  is 
driven  by  acoustic  motion.  Of  potentially  more  importance  to  propellant  combustion 
are  the  analogous  second-order  terms  which  appear  in  the  equation  for  the 
transport  of  thermal  energy  in  the  propellant  flame  zone.  While  these  thermal 
transport  effects  will  not  be  analyzed  in  this  study,  examples  of  their  importance  to 
effective  steady-state  heat  transfer  may  be  found  in  the  literature  (see,  for  example, 
Ref.  21).  The  effects  of  steady-state  enhancements  to  heat  transfer  from  acoustic 
streaming  effects  have  not  been  considered  with  respect  to  the  reactive  environment 
in  rocket  engines,  although  it  is  possible  that  these  effects  could  also  produce  a 
coupling  mechanism  for  instability. 

The  remainder  of  this  analysis  will  be  concerned  with  equation  (2)  and  the 
relative  importance  of  various  terms  within  the  acoust  ic  boundary  layer.  To  further 
simplify  this  equation,  it  will  be  assumed  that  the  density  in  the  region  of  interest  is 
approximately  constant  and  that  the  transverse  velocity,  v",  produced  by  propellant 
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response  is  negligible.  These  assumptions  would  necessarily  be  removed  in 
considering  a  more  comprehensive  propellant  velocity  response  analysis. 

The  magnitudes  of  the  terms  in  Equation  (2)  are  estimated  by  normalizing 
the  velocities  with  a0,  the  chamber  length  with  L,  the  chamber  radius  with  5  and 
utilizing  the  continuity  equation  and  the  assumption  of  a  standing  first  longitudinal 
acoustic  mode.  It  is  also  assumed  that  the  turbulence  shear  stress  may  be  scaled  by 
t"  =  ct  <p>  u"2,  where  I  ctl «  1 ,  The  respective  orders  of  the  terms  in  Equation  (2)  are 
then 


1,  <M>,  <M>,  <N1>,  (5,,'S)<M>.  <M>,  (a„/fSd)<Ms>,  1,  v/f5a-  cTa(>/(f5atturb)  M" }. 


The  order  unity  terms  (1)  and  (8)  recover  the  inviscid  acoustic  mode  solution 
in  the  central  region  of  the  duct.  A  considerable  simplification  results  if  <M>  is 
small,  since  terms  (2)  -  (6)  nur  be  dropped  without  fundamentally  altering 
boundary  layer  characteristics  retained  in  terms  (7),  (9)  and  (10).  (It  is  important  to 
note  that  the  maximum  response  to  the  surface  sensors  in  the  experiments  of 
Brown,  et  al.4  occurs  in  the  head  end  region,  where  <M>  is  <  0.1.)  Neglecting  for 
the  moment  the  transverse  convection  and  turbulence  terms,  it  is  seen  that  for  the 
viscous  term  to  be  of  order  unity  implies  the  Stokes  estimate, 


(3) 


and  is  the  smallest  possible  laminar  boundary-layer  height  of  the  problem.  The 
effects  of  injection  on  the  acoustic  boundary  layer  [provided  largely  by  term  (7)]  can 
be  substantial  and  even  dominate.  The  ratio  between  the  convection  and  viscous 
terms  in  Equation  (2)  is  of  the  order  10  under  conditions  of  interest. 


The  additional  shear  stress  provided  by  term  (10)  in  Equation  (2),  as  well  as 
the  analogous  additional  heat  flux  appearing  in  the  energy  equation,  result  from 
possible  turbulent  motions.  Further,  it  is  noted  that  term  (10)  is  the  only  retained 
term  which  depends  on  the  amplitude  of  the  acoustic  motion.  (The  acoustic 
amplitude,  M",  is  a  linear  scaling  parameter  of  all  other  terms  but  occurs 
quadratically  in  the  turbulent  shear  term),  thus  posing  the  question  of  a  possible 
relation  to  nonlinear  stability  phenomena. 
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As  noted  by  Merkli  and  Thomann12,  various  approaches  have  been  employed 
to  obtain  a  stability  or  transition  criterion  for  the  acoustic  boundary  layer.  In  the 
approach  adopted  here,  the  second-order  modeling  equations  employed  in  Ref.  5  are 
considered.  Preliminary  operations  and  assumptions  are  summarized  as  follows. 
Equations  for  the  velocity  correlations  u'  ju' j  are  contracted  to  obtain  an  equation 
for  the  turbulence  intensity  q2=u'iu’i-  Fluid  material  properties  are  assumed 
constant,  and  the  turbulence  (or  preturbulent  disturbance  level)  is  assumed  small  so 
that  third  and  higher  order  velocity  correlations  may  be  neglected.  It  is  also 
assumed  as  an  approximation  that  u’u'  =  auu<^.  v  v  =  awq2,  and  u  v  sa^q2, 
where  auu,  aw  and  auv  are  constants.  Substituting  these  relations  into  the  q2 
equation,  taking  the  mean  flow  to  be  negligible,  and  applying  boundary  layer 
assumptions  for  a  quasi-planar  condition  yields: 


au" 


p(  ^  +  u"  —  +  v"—  )  +  2p  q2(a  +  a  - 

p  at  ax  ay  uuay  uv  ay 


+  M- 


aV 

3y2 


(4) 


where  A  is  the  disturbance  macrolength  -  scale  and  A  is  a  constant  used  in  modeling 
low  turbulent  Reynolds  number  dissipation.  The  approximation  A  =  CA5a  (where  5a 
is  specified  by  Equation  (3)),  is  also  employed  to  evaluate  the  maximal  effective 
length  scale  appropriate  for  the  acoustic  boundary  layer. 


Equation  (4),  which  is  linear  in  q2,  may  be  order-of-magnitude  scaled  in  the 
same  manner  as  Equation  (2).  The  requirement  that  p(Dq2/Dt)  =  0  (for  neutral 
stability)  may  be  imposed  if  the  Reynolds  number  u  8a/u  <  0(1).  The  order  of 
magnitude  analysis  indicates  that  the  fifth  (production)  term  on  the  left  hand  side 
and  the  first  (dissipation)  term  on  the  right  hand  side  of  Equation  (4)  are  dominant, 
resulting  in  the  stability  criterion 


M"  -—A _ 

mcr  -  on.  a.. 


(-allv)C2,  a° 


(5) 
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Utilizing  the  turbulence  modeling  constants  specified  in  the  Appendix  of  Ref.  5  (i.e., 
A  =  3.25,  Ca  =  .17),  and  estimating  auv  =  -.15  from  fully  developed  flat-plate 
turbulent  boundary  layer  flows,  yields  an  estimate  of  K  =  750.  Merkli  and 
Thomann11  cite  prior  experimental  values  of  K  ranging  from  188  to  915,  and 
obtained  the  value  K  =  501  in  their  own  experiments.  The  constant  K  in  Equation 
(5)  is  related  to  the  Merkli  and  Thomann  constant  Ac  by  K  =  ACV (2ix)/2.  They 
speculated  that  the  variation  in  K  observed  in  prior  studies  could  be  caused  by 
variations  in  the  roughness  of  the  duct  surface,  and  demonstrated  that  disturbances 
caused  by  the  anemometer  probe  can  also  appreciably  affect  K. 

The  critical  Mach  number  given  by  Equation  (5)  is  shown  in  Figure  2  together 
with  the  data  of  Merkli  and  Thomann.  The  functional  dependence  is  correct,  and 
even  the  very  approximate  estimates  of  physical  constants  and  empirical 
parameters  yield  quantitative  ag;  >  ment  to  within  several  percent  of  the  data.  To 
the  authors'  knowledge,  the  derivation  of  an  approximate  transition  relation  from 
linearizing  this  type  of  complex  turbulence  model  is  novel. 
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0 

Figure  2.  Transition  to  turbulence  for  simple  acoustic  motion. 
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COMPUTATIONAL  ANALYSIS 


The  second-order  turbulence  closure  approach  developed  by  Donaldson  et 
al.22  and  Varma  et  al.23  was  implemented  in  Ref.  5  as  a  parabolized  model  for 
calculating  statistically  stationary,  compressible  transitional  flows  in  porous  walled 
ducts  with  large  injection  rates.  The  assumptions  concerning  the  gas  phase  are 
summarized  in  Appendix  A.  The  assumptions  and  order-of-magnitude  analysis  of 
the  previous  section  again  yield  a  parabolic  equation  system  for  the  present  acoustic 
boundary-layer  problem.  Computationally,  the  axial  convection  terms  pu  dg/dx 
(where  g  is  an  arbitrary  dependent  variable'  are  replaced  by  p  dg/dt  for  the  present 
case.  With  the  exception  of  dp/dx  in  the  momentum  equation,  all  other  axial 
derivatives  are  taken  to  be  null.  The  continuity  equation, 

—  +  1  —  [rV)(dv  +  p'v')3  =  0 
dt  rU  dr 

and  the  dp/dt  term  on  the  right  hand  side  of  the  energy  equation  are  retained, 
however,  for  subsequent  reactive  flow  investigations. 

The  parabolic  differential  equation  system  may  be  considered  in  the 
functional  form 


dt 


P~  +  pv 


r,  r)M 

dr 


=  J-i(r^g^)  +  Gg(f 

r"  dr  b  dr  b 


(6) 


where  g  =  i  u,h,  u  u'  ,  v’v' ,  w’w'  ,  uV  ,  h’u  ,  h'v'  ,  h  h  }T,  u,  v  and  w  are  the  axial, 
radial  and  circumferential  velocity  components,  v  is  the  index  for  planar  or 
axisymmetric  geometry,  p  the  density  and  h  =  cpT  is  the  specific  static  enthalpy. 
The  molecular  transport  coefficient,  p  represents  the  dynamic  viscosity  (g)  or  the 
thermal  conductivity  parameter,  a  =  k/cp,  as  appropriate  for  each  equation. 
However,  not  all  of  the  molecular  diffusion  terms  for  each  equation  in  the  system 
may  be  cast  in  the  form  shown  in  Equation  (2).  Those  that  do  not  conform  are 
implicitly  contained  within  the  complex  functions  0K,  which  also  represent  the 
sources,  cross-coupling  and  dissipation  terms  for  the  equations.  Although  the  more 
general  (  )  averaging  is  used  in  the  above  notation,  the  only  mean  velocity 

component  retained  is  v. 


The  pressure  along  the  duct,  p  =  vp>  +  p"(x,t),  is  specified  by  the  one¬ 
dimensional  standing  wave  solution 

p”  =  nma  <p>cos(nnx/L)cos(2  7tft)  (7) 

where  llma  is  the  maximum  relative  acoustic  pressure  amplitude,  and 
uma  =  IlmnVy.  Steep-fronted  waves  are  not  considered  in  this  investigation.  All 
nonreactive  flow  calculations  to  be  presented  were  performed  at  the  velocity 
antinode  of  the  first  longitudinal  mode  Material  properties  were  those  of  air  at 
standard  conditions,  or  at  elevated  pressures  where  noted. 

The  turbulence  length  scale  A,  in  this  model  is  algebraic  and  described^  by  a 
linear  variation  from  its  surface  value,  As  (As  =  0  for  smooth  walls),  to  a  plateau 
level  proportionate  (with  constant  Ca)  to  the  thickness  of  the  shear  flow.  For  the 
present  analysis  this  thickness  was  defined  as  Saa,  equal  to  the  height  above  the 
surface  where  the  boundary  asymptotes  to  99%  of  u"c.  For  injected  flows,  however, 
waves  are  convected  away  from  the  surface.  In  these  cases  the  first  zero-crossing  of 
u"  was  taken  as  the  effective  5aa.  The  final  values  of  5aa  used  in  the  length  scale 
expression  were  smoothed  by  integrating  in  time.  An  alternative  differential 
length-scale  equation  would  be  more  desirable,  assuming  that  problems  with  the 
low  turbulent  Reynolds  number  behavior  of  such  equations  could  be  circumvented. 

Boundary  conditions  at  the  duct  centerline  (or  eenterplane)  are  the  symmetry 
conditions,  which  for  the  posed  system  are 


—  [u,  h,  u  u'  ,  v'v  ,  w'w'  .  h'h'  ,  h'u'  ]  =  0  =u'v’  =  h'v'  . 
dr 

Due  to  the  final  form  of  the  governing  equations,  the  requirement  that  v  =  C  on  the 
centerline  (eenterplane)  is  necessarily  relaxed.  At  the  duct  surface,  the  static 
enthalpy  corresponding  to  a  given  temperature  is  specified,  and  the  mean  injection 
velocity  v3  is  prescribed.  The  no-slip  condition  implies  that  all  tangential  velocity 
components  and  their  correlations  are  zero.  Experiments  have  confirmed  that 
strong  injection  rates  through  a  porous  plate  can  produce  pseudo-turbulent 
disturbances  which  must  be  included  in  the  ensemble  average  of  the  v'v'  boundary 
condition.  The  parameter  o„  =  v'v'  fv  2  was  introduced  in  Ref.  5  to  account  for  this 
important  source  of  disturbance.  All  enthalpy  correlations  are  null  due  to  the 
prescribed  uniform  surface  enthalpy. 
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Gas-Phase  Combustion  Modeling 

The  gas  phase  equation  for  the  ensemble  mean  enthalpy,  under  the 
assumptions  outlined  in  the  previous  sections,  is  then 


Ip  v  +  p  v 


1  d 

rv  dr 


1  d  u  k  dh  \  dp 
ru  dr  c[<  dr:  dt 


+  M! 


2  i 


+  (A  +  B  Re^-U  AhR 

A“  i 


0) 


where  Ahg  to  is  the  energy  source  term  due  to  a  single-step  gas  reaction,  and  A  and 
B  are  constant  parameters  used  in  the  turbulence  model5-24.  For  the  present 
analysis  it  has  also  been  assumed  that  the  Lewis  number  is  unity,  and  that 
combustion  of  a  homogeneous  propellant  occurs  in  a  low  Mach  number  region.  The 
turbulent  enthalpy  correlations  above  the  flame  zone  are  necessarily  neglected  to 
avoid  contragradient  diffusion  which  can  result  in  numerical  instability.  The 
species  mass  fractions  are  expressed  in  terms  of  enthalpy,  due  to  Zeldovich 
similarity24,  even  though  the  cp/dt  term  in  the  enthalpy  equation  is  retained.  This 
assumption  is  based  on  an  order  of  magnitude  analysis  of  the  unsteady  mean 
enthalpy  equation,  valid  in  the  flame  zone,  which  yields  the  following  condition  for 
similarity: 


Pr(vs)x{rpj  ' 

I!  < 

f  R  kr 

The  above  expression  shows  that  the  relative  importance  of  the  dp/dt  term  in  the 
flame  zone  is  dependent  upon  the  pressure  ratio  of  the  <  sciPatmn  and  its  frequency 
in  addition  to  the  common  parameters,  us-uciated  with  quantifying  the 
characteristics  of  the  flame  zone.  Since  the  above  expression  is  satisfied  for  the 
example  calculations  to  be  discussed,  tie-  spe,  ms  mass  *V  ictions  needed  to  evaluate 


i 


needed  to  evaluate  the  reaction  rate  expression  are  expressible  in  terms  of  enthalpy 
using  the  similarity  approach  and  the  expression  for  w  is  then 


In  the  present  analysis,  o>  =  «  ('[\  j5)  since  it  has  been  shown  that  the  effects  of  near¬ 
surface  turbulent  reaction  rate  correlations  are  small  for  homogeneous 
propellants25. 


Unsteady  Heat  Conduction  in  the  Condensed  Phase 

The  unsteady  condensed-phase  heat  conduction  is  determined  through  a 
coupled  numerical  solution  of  the  gaseous  and  condensed  phases  employing  an 
iterative  procedure  based  on  the  gas/solid  interface  energy  balance.  The  condensed- 
phase  coordinate,  r\ ,  runs  from  zero  at  the  gas/solid  interface  to  -<»  in-  depth.  The 
gas-phase  coordinate,  y,  runs  from  zero  at  the  gas/solid  interface  to  +°°  on  the  scale 
of  the  thermal  boundary  layer.  The  assumptions  concerning  the  condensed  phase 
are  summarized  for  convenience  in  Appendix  A.  The  condensed  phase  heat 
conduction  equation,  in  the  absence  of  subsurface  reactions,  is: 


dT  .  dT  d2T 
—  +  r  - —  =  a„  — — 

dt  dr|  71  dr2 


(8) 


An  energy  balance  at  the  gas/solid  interface  yields  the  boundary  condition  at  r\  =  0: 


dTi  dTl 

■e'3yV 


(9) 


where  the  enthalpy  of  decomposition  is 
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t  o  =  y  ho  i, o 

-  Ld  <na'g  ~  nr 
a 

In  equation  (9),  the  assumption  of  large  chamber  radius  is  made  for  convenience. 
Equations  (8)  and  (9)  are  both  nonlinear  due  to  the  fact  that  the  surface  regression 
rate,  r  ,  is  a  strong  function  of  the  surface  temperature,  Ts.  The  regression 
(burning)  rate,  determined  by  condensed-phase  sublimation  or  evaporation,  is 
assumed  to  be  governed  by  an  Arrhenius  pyrolysis  law  of  the  form: 

m  =  |\  v„  =  pK  r  =  A  „  e  T'Vm  't«  (10) 

Finally,  the  in-depth  boundary  condition  is  simply: 

r\  -oo  ,  T  -> (ID 

where  Tn  is  the  constant  m-depth  propellant  temperature. 

Equations  (8)  through  (11)  completely  specify  the  unsteady  conduction  of  heat 
into  the  condensed-phase.  The  solution  is  obtained  by  performing  a  linear 
coordinate  transformation  on  equations  (8)  and  (9)  to  map  the  condensed  phase  onto 
the  computational  gas  phase  grid.  Equation  (8)  and  the  gas-phase  equations  (6)  are 
then  solved  numerically,  using  the  Crank  Nicolson  finite-difference  method,  and 
equation  (9)  is  used  to  iteratively  solve  for  the  interface  temperature.  The  linear 
transformation  from  rpspace  to  y-space  is: 

v  :  -  1  n  (12) 

Y 

where  =  krt  cp  /  (k(Ts)  cn)  is  a  scaling  parameter  used  to  approximately  match  the 
depth  of  the  combustion  wave  in  the  condensed-phase  to  the  flame  height  in  the  gas 
phase.  This  transformation  enables  the  condensed -phase  solution  to  take  full 
advantage  of  the  dynamically  adaptive  gas  phase  grid  and  ensures  consistent 
numerical  treatment.  Applying  the  transformation  to  equations  <H),  (9>,  and  (11  ) 
produces  the  following  system  of  fqmtinps 
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Numerical  Coupling  of  Condensed  and  Gas  Phase  Solutions 

The  gas/solid  interface  energy  balance,  Equation  (13b),  and  the  pyrolysis 
relation,  Equation  (10)  provide  the  coupling  between  the  condensed  and  gas  phase 
solutions.  The  coupled  solution  is  obtained  by  an  iteration  process  (here  referred  to 
as  a  substep)  at  the  new  time  step  as  follows: 


1)  Condensed  and  gas-phase  solutions  are  obtained  using  the  values  of 
the  gas/solid  interface  temperature  and  mean  injection  velocity 
specified  at  the  previous  temporal  step  or  substep. 


2)  An  updated  interface  mass  flux  is  specified  using  Equation  (13b)  is 
obtained  from: 
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where  n  and  n+1  are  the  nth  and  n+lst  substeps,  respectively. 


3)  equation  (10)  is  inverted  to  obtain  the  corresponding  surface 
temperature 
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and  the  updated  gas  phase  injection  velocity  is  determined  from 


(VS^  n+l  =!mS'!!n+l  -''(Pgls 

the  condensed  and  gas  phases  are  solved  once  again,  using  the  updated 
values,  and  the  process  is  repeated  until  convergence  is  achieved 
(typically  in  two  or  three  iterations). 


STABILITY  OF  ISOTHERMAL  ACOUSTIC  BOUNDARY  LAYERS 


The  initial  profile  for  the  axial  velocity  component  was  taken  to  be  the 
analytical  solution  (at  t  =  0)  of  Sex’  and  tJchida  fin  Schlichting25)  for  oscillatory  flow 
in  a  cylindrical  duct.  Initial  profiles  of '.he  normal  Reynolds  stresses  were  taken  to 
be  isotropic,  and  proportional  to  the  square  of  the  initial  velocity  profile.  The  peak 
initial  disturbance  level  qrnax*1  =  0,yVuma,  was  assumed  to  be  0.01  for  all 
calculations  reported  here. 

The  implicit  numerical  procedure  and  adaptive  grid  are  described  in  prior 
studies^.  Except  where  nou  d,  approximately  50  time  steps  per  period  were 
employed.  The  number  of  radial  spatial  nodes  within  the  boundary  layer  varied 
from  about  40  for  laminar  flows  to  about  75  for  turbulent  cases. 

Normalized  profiles  of  the  axia1  city  ;r>  the  laminar  acoustic  boundary 
layer  (f  =  1000  Hz,  5a  =  180  urrf  aie  shown  for  various  times  within  a  period  in 
Figure  3.  The  results  display  the  classical  phase  shift  and  Richardson  annular 
effects  due  to  viscosi;  r.  These  calculations  were  performed  with  100  time  steps  per 
period  and  produce  a  relatively  small  computational  error,  principally  outside  the 
boundary  layer  and  near  the  maxima  in  the  acoustic  velocity.  However,  since  a 
large  number  of  calculations  were  needed  for  the  transition  studies,  the  number  of 
time  steps  per  period  was  halved  for  those  calculations,  and  the  error  in  the 
calculated  acoustic  velocity  was  then  approximately  doubled. 

Figure  4  shows  the  axial  velocity  profiles  for  a  fully  (cyclically)  turbulent 
boundary  layer  at  100  Hz  with  llma  -  .125.  There  is  a  noticeable  diffusion  of  the 
Richardson  effect  and  a  pronounced  increase  in  the  velocity  gradient  near  the 
surface.  However,  a  protracted  law-of-thc -wall  (logai  ithmic)  region  is  not  observed. 
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Figure  3.  Laminar  velocity  profiles  in  the  acoustic  boundary  layer 
without  injection. 

Figure  5  shows  normalized  turbulence  intensity  profiles,  q/uma,  in  the 
boundary  layer  for  the  same  conditions  stated  for  Figure  4.  The  calculated  peaks  in 
turbulence  level  are  comparable  to  those  calculated  for  steady  state  flows  (=  15%). 
To  within  computational  accuracy,  these  peaks  are  symmetric  with  respect  to  ±rcft. 
A  slight  phase  lag,  on  the  order  of  5-10  degrees,  exists  between  the  intensity 
maxima  and  the  rectified  acoustic  velocity  I  u”c  I .  This  lag  increases  at  higher 
frequencies.  Also  note  that  the  intensity  decreases  substantially,  but  is  not 
predicted  to  vanish  at  the  zero-crossings  of  u  ”c. 

Figure  6  shows  turbulent  velocity  profiles  for  a  Stokes  layer  with  injection  at 
f  =  1  00  Hz,  II  =  0.080,  and  mean  injection  velocity  <v>  --  1  m/s.  A  pronounced 
Richardson  effect  is  evident,  surpassing  even  the  laminar  noninjected  results.  The 
estimated  value  of  the  effective  boundary  layer  thickness  for  this  case  is  5.5  mm 
(approximately  y/6a  =  14  in  the  figure).  The  convective  wave-train  behavior  shown 
is  also  evident  in  the  laminar  injected  Stokes  layer  calculations  of  Click  and  RenieS 
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Figure  4.  Turbulent  velocity  profiles  in  the  acoustic  boundary  layer 
without  injection. 
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Figure  5.  Turbulence  intensity  profiles  in  the  acoustic  boundary  layer 
without  injection. 


19 


O  25  50  75  1  OO 


Normalized  Distance  From  Surface,  y/'<5 

a 

Figure  6.  Turbulent  velocity  profiles  in  the  acoustic  boundary  layer 
with  injection. 
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Figure  7.  Turbulence  intensity  profiles  in  the  acoustic  boundary  layer 
with  injection. 


Figure  7  shows  the  corresponding  normalized  turbulence  intensity  profiles. 
Due,  in  part,  to  the  enhanced  Richardson  effect,  very  large  intensity  levels  (=  30%) 
are  predicted.  As  turbulence  is  produced  in  the  inner  layer  at  a  later  time  in  the 
cycle,  the  velocity  wave  formed  from  an  earlier  cycle  produces  a  harmonic  wave  in 
the  intensity  profile  (ft  =  5.125).  A  long  tail  is  predicted  due  to  large  decay  times 
and  the  fact  that  turbulence  is  beginning  to  accumulate  at  the  computational 
centerline  (8  =  5  cm)  where  the  symmetry  condition  is  employed.  In  an  actual  duct 
flow  environment,  the  mean  axial  flow  would  convect  these  tails  downstream  as 
they  interact  with  mean  flow  induced  turbulence. 

Calculated  results  for  the  critical  acoustic  Mach  number  at  transition, 
M  "Cr  =  ncr/y,  are  shown  in  Figure  2  together  with  the  approximate  relation  and  data 
for  non-injected  Stokes  layers.  The  functional  dependence  is  similar,  but  the 
proximity  of  the  calculations  to  the  K  =  750  line  is  coincidental.  Note,  for  example, 
that  although  not  shown  here,  the  numerical  results  display  a  dependence  on  the 
initial  amplitude  of  turbulence  assumed.  “Transition”  in  the  calculated  results  was 
determined  by  monitoring  the  growth  or  decay  of  initial  turbulence  from  the  initial 
value  over  several  cycles.  For  the  cases  with  injection,  the  growth  rates  were  not  as 
strongly  dependent  on  as  were  those  for  the  non-injected  cases. 

Critical  Mach  numbers  for  transition  as  a  function  of  the  normalized  mean 
injection  velocity,  <v>/(fv)1''2,  are  shown  in  Figure  8  for  three  frequencies,  f  =  100, 
300  and  1000  Hz.  A  pronounced  decrease  in  M"cr  is  predicted  with  increasing 
injection  velocity,  so  that  all  curves  exhibit  a  minimum  that  is  frequency  dependent. 
The  minima  are  almost  a  factor  of  three  below  the  M”cr  for  noninjected  Stokes 
layers.  As  an  example,  the  100  Hz  non-injected  Stokes  layer  is  calculated  to  have 
M"cr  =  .083.  For  the  injected  Stokes  layer,  the  minimum  value  ofM"cr  =  .025  occurs 
at  <v>  =  0.25  m/s,  while  M'  cr  =  .048  at  <v>  =  1.0  m/s  (the  termination  of  the  line  at 
this  frequency). 


21 


Figure  8.  Effect  of  injection  on  the  stability  of  the  acoustic  boundary 
layer. 

Even  more  informative  are  the  same  results  for  transition  normalized  not  as 
a  critical  acoustic  Mach  number,  but  as  a  critical  acoustic  Reynolds  number, 
UmaAfu)1^2.  In  this  case,  all  three  curves  for  frequency  in  Figure  8  collapse  onto  a 
single  line  (to  within  computational  precision),  shown  as  the  upper  line  in  Figure  9. 
This  Figure  also  shows  results  for  finite  levels  of  the  surface  disturbance  parameter 
ov  =  0.035  and  0.07  for  two  different  transition  criteria  to  be  discussed.  The 
accompanying  value  of  As  assumed  in  the  calculations  was  3x1 0"^  m.  These  values 
are  viewed  as  realistic  and  perhaps  conservative.  It  is  not  surprising  that  this  finite 
disturbance  effect,  which  continuously  “feeds”  the  injected  layer,  is  predicted  to 
yield  further  reduction  in  critical  acoustic  Reynolds  number,  analogously  to  quasi¬ 
steady  flows. ^  As  a  specific  example,  the  minimum  critical  Mach  number  for  the 
100  Hz  case  at  <v>  -  .20  m/s  is  decreased  to  approximately  one  percent  for  ov  = 
0.035. 
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Figure  9.  Effects  of  injection  and  injection  disturbance  level  on  the 
stability  of  the  acoustic  boundary  layer. 

As  mentioned  earlier,  it  was  more  difficult  to  determine  the  value  of  M"cr  for 
cases  with  injection  and  finite  levels  of  av  present.  The  first  criterion  employed 
consisted  of  monitoring  the  maximum  turbulence  level  qmax^ma  occurring  at  any 
height  within  the  boundary  layer.  If  this  maximum  attained  a  value  of  0.07  or 
greater  within  3  cycles,  the  flow  was  considered  turbulent.  Using  this  criterion,  the 
ov  =  0  lines  were  computed  along  with  the  alternating-dashed  lines  for  ov  =  0.035 
and  0.07.  While  this  criterion  was  satisfactory  for  av  =  0  cases,  it  proved  to  be 
unsatisfactory  for  cases  with  finite  levels  of  av. 

To  illustrate  this  effect,  Figure  10  shows  the  variation  of  turbulence  intensity 
vs  normalized  time  at  y  =  5a,  for  three  values  of  n,  <v>  =  0.25  m/s  and  cv  =  0.  It  is 
evident  from  the  curves  in  Figs.  10a  and  10b  that  a  very  small  increase  in  n  is  able 
to  cause  the  flow  to  transition  to  a  cyclically  turbulent  condition. 


Figure  10  (a-c).  Time  variation  of  turbulence  development  for  the  injected 
acoustic  boundary  layer  (av  =  0,  y  =  5a). 

When  finite  levels  of  ov  are  included,  the  determination  of  ncr  becomes  more 
complicated.  Figures  11a,  b,  and  c  are  results  for  <v>  -  0.25  riv's  and  cv  =  0.035. 
Figure  11a  is  a  case  for  which  the  turbulence  bursts  to  a  q/Uma  level  of 
approximately  0.24  but  then  dissipates  over  10  cycles  until  only  the  low  level 
turbulence  due  to  ov  is  still  present.  Figures  1  lb  and  c  illustrate  the  increase  in  II 
required  to  achieve  the  cyclic  behavior  which  we  considered  to  be  characteristic  of  a 
turbulent  case.  This  criterion  for  tiensition  was  used  to  calculate  the  dashed  curves 
for  and  0.07  in  Figure  9. 

To  recapitulate  the  results  of  the  two  transition  criteria,  the  alternating 
dashed  lines  in  Figure  9  qualitatively  represent  a  “burst”  or  transient  criterion  for 
transition,  while  the  dashed  lines  represent  the  particular  'long-time)  criterion.  For 
the  latter  case,  regions  to  the  right  of  the  nearly  vertical  line  segments  are  indicated 
to  be  turbulent.  The  difference  between  the  two  criteria  is  in  part  explained  by 
Equation  t4),  which  must  be  augmented  by  the  term  p  '"v>  when  is  finite. 
This  can  result  in  three  different  ( urhiPencc  time  -calcs  being  competitive  under 
some  conditions. 
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Figure  11  (a-c).  Time  variation  of  turbulence  development  for  the  injected 
acoustic  boundary  layer  (ov  =  0.035,  y  =  5a). 


REACTIVE  ACOUSTIC  BOUNDARY  LAYER  -  PRESSURE  NODE 

The  thermochemical  properties  and  reference  operating  conditions  of  the 
baseline  solid  propellant,  presented  in  Appendix  B,  were  chosen  to  represent  the 
nominal  properties  of  homogeneous  propellants  and  were  also  used  in  the  erosive- 
burning  study  of  Ref.  26.  Calculations  presented  were  performed  using  the  baseline 
propellant,  which  has  an  exothermic  surface  reaction  (L°  =  -4.184  x  105  J/kg),  and 
an  identical  propellant  with  an  endothermic  surface  reaction  (L°  =  837  J/kg).  To 
reproduce  the  reference  operating  conditions,  the  pre-exponential  constant  in  the 
gas  phase  reaction  rate  was  calculated  to  be  8.2x1 09  kg^m^-’s^K8^  for  the 
endothermic  propellant.  Steady-state  calculations  performed  over  a  broad  range  of 
pressures  yielded  a  normal  rate  of  regression  r  =  r  (p/p*f  with  r  and  n  as  indicated 
in  the  Table.  Transition  calculations  were  performed  using  50  time  steps  per  cycle 
and  two  iterations  per  step  while  propellant  response  calculations  were  performed 
using  100  steps  per  cycle  and  throe  iterations  per  step.  Unless  stated  otherwise,  all 
propellant  responses  were  obtained  at  an  acoustic  frequency  of  1000  Hz  and  five 
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percent  acoustic  pressure  ratio,  IIma  =  0.05.  The  disturbance  level  of  the  gas 
transpiring  from  the  propellant  surface,  ctv.  was  assumed  to  be  zero  in  all 
calculations. 

Shown  in  Figure  12  are  the  normalized  steady-state  temperature  profiles 
obtained  in  the  gas  and  condensed  phases  for  the  endothermic  and  exothermic 
propellant  formulations  at  reference  conditions,  normalized  by  the  approximate 
theoretical  gas  phase  flame  height,  8f,  (note  that  the  approximate  8f  scaling  used 
here  is  independent  of  L").  The  results  indicate  that  the  exothermic  formulation 
provides  an  anticipated  increased  gas  phase  flame  height,  increasing  the  possibility 
of  acoustically  induced  turbulence  interacting  with  the  flame  zone.  The  coincident 
condensed  phase  profiles  are  shown  plotted  in  the  gas  phase  computational 
coordinates  and  are  therefore  scaled  by  ^  =  4.0  in  all  calculations.  Thus,  ;f  the 
scaling  were  not  employed,  the  condensed  phase  combustion  wave  would  be 
approximately  four  times  the  length  of  the  flame  zone  and  the  condensed  phase 
solution  would  not  efficiently  use  the  dense  grid  (generally  55  points)  provided  by 
the  dynamically  adaptive  grid  within  the  combustion  zone. 


Steady  state  t<  mperniure  profiles. 


Figure  12. 


T'igure  13  presents  the  normalized  temperature  and  turbulent  velocity 
profile  vs.  distance  from  the  burning  surface.  The  abscissa  is  normalized  by 
5a  =  10  jtm.  The  profiles  a>-o  given  °or  various  times  within  the  acoustic  cycle  for 
the  ba:  ?line  propellant  with  an  equi  /alent  roughness  of  ks  =  150  pm.  The  abscissa 
is  pres  r.ted  as  a  logarithmic  scale  to  reveal  the  details  in  the  near  surface  flame 
regioa  The  results  obtained  show  the  phase  shift  and  Richardson  annular  effect 
due  to  viscosity  and  injection  that  were  obtained  in  the  previous  section  for 
isothermal  Stokes  layers.  The  acoustic  velocity  is  negligible  in  the  region  of  the 
flame  zone  at  this  intermediate  frequency.  Figure  14  shows  the  corresponding 
normalized  temperature  profilet,  anti  turbulence  intensity  profiles.  The  peaks  in 
turbulence  intensity  within  the  cycle  correspond  approximately  to  the  large  velocity 
gradients  produced  by  the  Richardson  annular  effect.  The  results  indicate  that  the 
turbulence  intensity  obtains  magnitudes  of  one  to  two  percent  within  the 
combustion  region  at  certain  points  within  the  acoustic  cycle,  and  results  in 
propellant  combustion  responds-  to  enhanced  heating  rates  due  to  turbulent 
diffusion. 

In  a  previor  ,  .cion,  it  was  found  that  for  a  smooth  surface  (ks  =  0)  without 
surface  injection  disturbances  iov  =  0),  an  acoustic  transition  locus  could  be 
expressed  m  terms  of  axial  acoustic  velocity  Reynolds  number  as  a  function  of  an 
injection  velocity  acoustic  Reynolds  number.  Figure  15  shows  the  previously 
reported  results.  The  critical  amplitude  locus  is  quite  sensitive  to  the  method  used 
to  introduce  “seed”  turbulence  into  the  boundary  layer  and  to  other  conditions  such 
as  the  relative  height  within  the  boundary  layer  at  which  amplification  is 
determined.  Figure  15  shows  additional  calculations  for  an  isothermal  boundary- 
layer  which  utilize  an  alternative  means  of  introducing  the  initial  disturbance.  (The 
difference  between  the  two  methods  is  that  the  earlier  results  introduced  a  lc/r  peak 
disturbance  level  within  the  Stokes  layer  thickness  5a,  while  the  present  results 
introduce  the  same  19f  disturbance  level  over  the  much  larger  actual  injected 
boundary  layer  thickness,  <5aa). 
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Figure  13.  Comparison  of  temperature  and  turbulent  velocity  profiles  at 
various  times  within  a  cycle. 
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Figure  14.  Comparison  of  corresponding  temperature  and  turbulence 
intensity  profiles  at  various  times  within  a  cycle. 
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Figure  15.  Critical  acoustic  Reynolds  number  required  for  transition. 

Using  this  new  and  perhaps  more  realistic  method,  transition  calculations 
were  also  performed  for  the  baseline  propellant  at  frequencies  of  100  and  1000  Hz 
over  p’-ossures  ranging  from  10  to  100  atmospheres.  The  transition  points  obtained 
for  propellant  combustion  are  in  good  agreement  with  the  isothermal  calculations. 
Thus,  the  appropriate  scaling  involves  flame  or  edge  quantities  only,  and  indicates 
that  the  acoustic  flowfield  (for  a  smooth  surface)  can  be  approximated  as  an 
isothermal  flow  at  the  steady-state  flame  temperature  in  determining  transition. 

Figure  16  shows  the  calculations  presented  in  Fig  15,  for  the  endothermic 
propellant  along  with  calculations  for  the  exothermic  propellant  at  an  acoustic 
frequency  of  100  Hz,  plotted  as  a  function  of  the  critical  acoustic  pressure  ratio,  ncr. 
vs.  the  mean  pressure  normalized  by  p*.  The  critical  pressure  ratio  is  found  to  be 
insensitive  to  the  heat  released  from  the  surface  reaction,  although  the  actual 
response  of  the  propellant  is  dependent  on  L”.  This  results  from  the  fact  that,  for  a 
smooth  wall,  transition  is  initiated  above  the  thin  combustion  zone.  The  critical 
pressure  ratio  is  found  to  decrease  with  increasing  mean  pressure  in  all  cases, 
indicating  that  a  propellant  which  may  be  stable  with  respect  to  a  specific  acoustic 
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node  at  a  known  mean  pressure  may  encounter  acoustic  transition  at  that  mode 
/hen  operating  at  an  elevated  mean  pressure.  This  phenomena  may  explain 
tability  problems  encountered  in  the  development  of  solid  rocket  motors  with 
mconventionally  high  chamber  pressures. 
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figure  16.  Variation  of  critical  pressure  ratio  required  for  transition  as  a 
function  of  mean  chamber  pressure. 

Figure  17  shows  the  response  of  the  baseline  endothermic  propellant  to  the 
presence  of  the  turbulent  acoustic  boundary-layer  for  various  surface  roughnesses. 
The  normalized  propellant  mass  flux  'burn  ram;  ;<  plotted  vs  normalized  time,  t.  A 
liagram  of  the  variation  of  the  acoustic  edge  velocity  and  pressure  is  included  at  the 
op  of  the  figure  The  pressure  vaiiad'  n  is  cm  »<  ■•■!  fm-  an  avia!  location  in  the  head- 
>nd  region  of  a  chamber  In  the  att-end  region,  tin.*  pressure  would  be  shifted  180 
legrees  in  phase  with  respect  to  the  pressure  curve  shown.  The  case  of  a  smooth 
vail,  ks=  0  m,  required  a  10  percent  pre—  ere  ratio  to  obtain  the  response  shown. 
The  cases  in  which  a  finite  value  of  surf-m."  roughness  was  prescribed  were 
:alculated  using  a  five  percent  pi*  ssure  i  atio  The  smooth  wall  produced  only  a 
slight  response  in  the  propellant  !>u*-n  ecr  !liat  >•-'  i  -  >  > fi >  ol  with  respect  to  the 
velocity  field  and  which  is  negligible  t.vc>  a  larg*  pori  im>  of  each  half  cycle.  The 
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calculations  in  which  finite  levels  of  surface  roughness  are  included  reveal  several 
interesting  features.  All  cases  exhibit  the  rectified  response  along  with  a  time-mean 
augmentation  (“D-C  shiii”)  of  the  propellant  burn  rate,  about  eight  percent  for 
ks  =  150pm,  five  percent  for  125pm,  and  two  percent  for  the  100pm  surface 
roughness.  The  relative  phase  of  the  response,  with  respect  to  the  acoustic  velocity 
and  pressure  field,  varies  with  surface  roughness  and  may  provide  a  mechanism  by 
which  the  acoustical  energy  can  he  reinforced  or  damped. 

Figure  18  compares  the  e fleet  surface  decomposition  energy  on  the 
propellant  response.  The  calculations  were  performed  for  the  endothermic  and 
exothermic  propellant  formulations  with  -  150pm,  and  IT  —  0.05,  in  both  cases. 
The  exothermic  propellant  has  a  mean  shift  in  burning  rate  that  is  approximately 
three  times  as  large  as  that  occurring  in  the  endothermic  case. 

The  propellant's  respon  i  to  the  acoustic  boundary  layer  is  further 
investigated  using  a  conventional  Fast  Founer  Transform  (FFT)  analysis  over 
several  cycles.  The  cases  examined  are  identical  to  the  calculations  with  finite 
roughness  presented  in  Fig  17  and  discussed  above.  Figures  19  and  20  present  the 
results  of  the  analysis  in  terms  of  the  relative  amplitude  and  frequency  content  of 
the  response  and  the  relative  phase  of  the  responses  at  each  frequency, 
respectively.  The  large  dashed  curve  is  the  driving  frequency  of  the  external 
acoustic  velocity,  the  small  dashed  curve  is  the  response  of  the  propellant  with  a 
smooth  wall  (the  response  has  been  graphically  shifted  +300  Hz  to  facilitate  the 
comparison),  and  the  solid  curve  is  the  response  with  ks  =  150  pm.  The  frequency 
content  and  relative  amplitudes  of  the  response  with  ks  =  150  pm  is  representative 
of  the  cases  with  finite  surface  roughness.  The  results  presented  in  Fig  19  have,  in 
each  case,  been  normalized  by  the  largest  amplitude  obtained  for  that  calculation 
after  removal  of  the  zero-frequency  component.  It  can  be  seen  that  in  both 
calculations  the  primary  response  occurs  at  the  second  harmonic  with  significant 
harmonic  content  occurring  at  the  lower,  even  harmonics.  Only  even  harmonics  are 
obtained  for  this  case  since  axial  mean  flowr  is  neglected  and  calculations  are 
performed  at  a  pressure  node.  The  relative  phase  of  the  propellant  response, 
presented  in  Figure  20,  is  seen  to  vary  substantially  from  harmonic  to  harmonic.  It 
can  further  be  seen  that  at  the  second  harmonic,  2000  Hz,  the  phase  lag  of  the 
response  increases  with  increasing  surface  roughness  (excluding  the  smooth  wall 
case  due  to  the  large  pressure  ratio  required  to  obtain  the  response). 
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Figure  17. 


Figure  18. 


Normalized  Time,  r  (cycles) 


Variation  of  propellant  response  as  a  function  of  time  for 
various  surface  roughness  heights. 
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Variation  of  propellant  response  as  a  function  of  time  for 
endo-  and  exothermic  surface  decomposition  energies. 


Figure  19. 


Figure  20. 
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Harmonic  content  of  the  propellant  responses  shown  in 
Figure  16  for  a  smooth  and  rough  wall  (ks  =  0  and  k8  =  150  gm, 
respectively) 


Relative  phase  angle  of  the  responses  shown  in  Figure  16. 
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COMBINED  RESPONSE  AT  VARIOUS  AXIAL  LOCATIONS 


The  thermophysical  properties  and  reference  operating  conditions  of  the  solid 
propellant  employed  in  the  previous  section  were  also  used  to  investigate  the 
combined  response  and  acoustic  erosion  at  various  axial  locations.  All  calculations 
presented  were  performed  using  100  time  steps  per  cycle  and  three  iterations  per 
time  step.  When  calculating  the  acoustic  erosion  of  the  propellant  at  various  axial 
locations,  the  integration  of  the  propellant  response  was  performed  in  the  7th  cycle 
of  oscillation  after  the  introduction  of  the  seed  turbulence  in  order  to  let  the 
propellant  combustion  reach  a  cyclically  steady-state.  The  integration  required  to 
calculate  the  acoustic  erosion  due  to  acoustic  velocity  oscillations  alone  were 
performed  in  the  3rd  cycle  of  oscillation  after  the  seed  turbulence  was  first 
introduced  since  the  flow  would  reach  the  cyclically  steady-state  condition  in 
approximately  two  cycles.  All  calculations  presented  were  performed  at  an  acoustic 
frequency  of  1000  Hz  and  five  percent  acoustic  pressure  ratio,  nn,a  =  0.05,  unless 
otherwise  stated.  The  calculated  acoustic  boundary  layer  heights  were  5aa  =  1mm 
and  0.6mm  for  an  acoustic  frequency  of  f=1000  Hz  and  2000  Hz,  respectively. 

The  condensed-phase  solution  was  calculated  on  the  gas-phase  computational 
grid  after  being  scaled  by  ^  =  4.0  to  use  the  dense  grid  provided  by  the  combustion 
zone.  The  flame  zone,  which  has  a  height  of  approximately  10pm,  generally  had  35 
grid  points  within  it  while  the  acoustic  boundary-layer  and  remainder  of  the 
computational  domain,  which  extends  to  a  height  of  5cm,  contained  approximately 
60  and  180  grid  points,  respectively.  The  disturbance  level  of  the  gas  transpiring 
from  the  propellant  surface,  ctv,  was  assumed  to  be  zero  in  all  calculations.  The 
acoustic  erosive  burning  is  calculated  by  integrating  the  instantaneous  propellant 
regression  rate,  r,  over  one  cycle  of  oscillation  to  obtain  the  long  time  mean  of  the 
regression  rate,  <-r  >,  which  is  then  normalized  by  the  normal  burning  rate  of  the 
propellant,  rn 

Shown  in  Figure  21  are  the  steady-state  temperature  profiles  in  the  gas  and 
condensed-phases  normalized  by  the  steady-state  flame  temperature,  Tf,  as  a 
function  of  distance  from  the  burning  surface  normalized  by  the  approximate 
theoretical  laminar  flame  height.  Of.  The  condensed-phase  temperature  profile  is 
plotted  in  the  gas-phase  computational  coordinates  and  has  therefore  been  scaled  by 
=  4.0  as  it  is  in  all  calculations  The  calculated  profiles  are  m  good  agreement 


with  the  approximate  theoretical  flame  height  and  theoretical  condensed-phase 
thermal  boundary-layer  depth. 


Distance  From  Burning  Surface,  y/<5f 

Figure  21.  Normalized  steady-state  temperature  profiles. 

Figure  22  is  a  comparison  of  the  gas-phase  temperature  perturbations  due  to 
acoustic  pressure  oscillations  (laminar  case)  calculated  both  with  and  without  the 
assumption  of  Zeldovich  similarity.  The  latter  case  being  equivalent  to  the  equation 
system  used  in  Tien's  analysis.27  The  calculations  were  performed  at  x/L=0.375 
with  nma=0.05  (which  yields  a  locrl  acoustic  pressure  ratio  of  1.913  percent 
(ri|=0.01913)  and  are  found  to  be  virtually  identical.  The  perturbations  presented 
are  defined,  following  the  analysis  of  Ref.  27,  as  the  difference  between  the 
instantaneous  value  and  the  steady-state  value  of  the  variable  normalized  by  a 
reference  condition  (normally  a  value  occurring  at  the  edge  of  the  flame  zone)  and 
the  local  acoustic  pressure  ratio.  Various  profiles  throughout  one  cycle  of  oscillation 
are  plotted  as  a  function  of  normalized  distance  from  the  burning  surface  to  a 
distance  of  1000  theoretical  flame  heights.  The  perturbations  qualitatively  agree 
with  the  results  of  Ref.  27  and  verify  the  oscillatory  behavior  of  the  wave  train 
exterior  to  the  flame  zone.  A  direct  comparison  with  the  results  of  Ref.  27  was  not 
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possible  since  thermochemical  parameters  of  a  representative  solid  propellant  were 
selected  for  this  study,  where  as  the  nondimensional  parameters  of  Ref.  27  produce 
inordinately  low  regression  rates  and/or  high  frequencies.  This  would  imply  that  the 
(dimensional)  frequencies  required  for  a  direct  comparison  with  the  results  of  Ref. 
27  be  the  order  of  106  Hz.  At  frequencies  of  this  magnitude,  many  of  the 
assumptions  of  both  analyses  would  be  invalid.  Calculations  up  to  104  Hz  were 
performed  which  appeared  to  be  in  qualitative  agreement  with  the  perturbation 
results  of  Ref.  27. 
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Figure  22.  Comparison  of  normalized  temperature  perturbations  with 
and  without  the  similarity  assumption. 


The  curves  presented  in  Figure  23  are  nf  the  laminar  and  turbulent  gas- 
phase  temperature  perturbations  occurring  at  various  points  within  one  cycle  of 
oscillation  plotted  as  a  function  of  normalized  distance  from  the  burning  surface 
corresponding  to  the  gas  phase  flame  /.one  plotted  in  Figure  21  The  calculations 
were  performed  at  x/L-0.375,  with  I!|=0  itho  value  of  1 1  at  x/L)  and  f=l 000 

Hz.  The  low-amplitude,  lamina'-  flow  prrtuH  at  inn-  oscillate  symmetricallv  about 
the  steady  state  solution,  as  expected  '  iwrlinnul  periurhation?  are  found  to  all 


be  positive,  accounting  for  the  increase  of  approximately  5  percent  in  the  mean 
regression  rate  of  the  condensed-phase  (DC-shift)  at  this  location,  and  are 
substantially  larger  in  magnitude  than  the  corresponding  laminar  perturbations.  It 
is  noted  that  the  “turbulent”  perturbations  presented  are  actually  the  combined 
response  due  to  both  turbulent  flow  and  harmonic  pressure  fluctuations,  but  will  be 
referred  to  as  the  “turbulent”  perturbations  in  subsequent  text  for  brevity. 


Distance  Drum  Burning  Surface,  y/c5^ 

Figure  23.  Normalized  laminar  and  turbulent  gas-phase  temperature 
perturbations. 

Figure  24  compares  the  laminar  and  turbulent  gas  phase  normal  velocity 
perturbations  directly  corresponding  to  the  temperature  perturbations  presented  in 
Figure  23.  The  laminar  perturbations  have  been  scaled  in  magnitude  by  a  factor  of 
5  for  presentation.  It  is  evident  that  the  pressure  response  is  dominated  by  the 
effects  of  turbulence  produced  by  the  acoustic  motion.  The  turbulent  normal 
velocity  perturbations  at  the  surface  and  edge  of  the  flame  zone  are  found  to 
oscillate  about  a  mean  shift  corresponding  to  the  shift  in  the  regression  rate  of  the 
condensed-phase. 


The  laminar  and  turbulent  perturbations  of  the  condensed-phase 
temperature,  directly  corresponding  to  the  curves  presented  in  Figs.  21 , 23  and  24, 
and  plotted  in  the  gas-phase  computational  space  are  shown  in  Figure  25.  The 
magnitude  of  the  laminar  perturbations  have  been  -h  a  led  by  a  factor  of  5  for 
presentation.  The  laminar  perturbations  do  not  oscillate  about  the  steady-state 
profile,  as  the  laminar  gas  phase  perturbations  do,  but  oscillate  about  an 
asymmetrical  mean  perturbation.  This  asymmetry  arises  from  the  variation  of  the 
regression  rate  land  therefore,  convection  i  of  the  condensed -phase  throughout  the 
cycle  of  oscillation  in  combination  with  the  oscillator',  nature  of  the  temperature  at 
the  gas-solid  interface  The  turbulent  perUmbations  are  found  to  be  on  the  order  of 
ten  times  larger  than  the  corresponding  laminar  perturbations.  The  extremely 
asymmetrical  behavior  of  the  turbulent  profiles  is  due  primarily  to  the  increase  in 
the  mean  regression  rate  of  the  condensed-phase  at  this  location,  which  tends  to 
reduce  the  depth  of  the  ’hernial  hounda?-y  lay-r.  The  laminar  and  turbulent 
acoustic  velocity  profiles  and  the  corresponding  turbulence  intensity  profiles  are 
essentially  identic  al  to  'he  a  p.-u  dm,'  fw  ■ .  i-  t h  p ■ .  \  mu-  section. 


Figure  25.  Normalized  laminar  and  turbulent  condensed-phase  temper¬ 
ature  perturbations. 

Results  have  been  obtained  for  the  response  of  the  propellant  regression  rate 
to  conditions  which  would  be  representative  of  those  existing  in  a  center-vented  T- 
burner,  or  in  the  head-end  region  of  a  solid  rocket  motor  where  the  mean  axial  Mach 
number  is  low,  when  either  are  experiencing  acoustic  oscillations.  For  the 
fundamental  mode,  the  acoustic  environment  varies  from  the  pure  pressure 
oscillation  at  either  end  to  the  pure  velocity  oscillations  in  the  center  of  the  chamber 
at  the  pressure  node,  and  a  combination  of  both  at  any  point  between.  Only  results 
for  the  portion  of  the  duct  from  x/L=0  to  x/L=0.5  will  be  presented  since  in  the 
absence  of  a  mean  axial  flow  the  response  will  be  cyclically  identical  in  either  end  of 
the  duct. 

The  combined  response  of  the  condensed  phase  regression  rate  at  various 
positions  along  the  chamber  are  shown  in  Figure  26.  The  top  portion  of  the  plot 
contains  a  diagram  of  the  corresponding  acoustic  pressure  oscillation  and  the  two 
possible  external  acoustic  velocity  oscillations  which  could  exist  at  each  location 
along  the  chamber  (the  diagram  shows  only  the  relative  phases  of,  and  is  not 
intended  to  show  the  amplitude  relationship  between,  the  acoustic  velocity  and 
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pressure).  The  calculations  were  performed  for  Mma  =  0.05  and  ks  =  150pm.  The 
response  of  the  condensed-phase  repression  rate  at  the  end  of  the  chamber,  x/L=0,  is 
sinusoidal  and  lags  the  pressure  by  approximately  9  degrees.  At  a  quarter  point  in 
the  duct,  x/L=0.25,  the  turbulent  acoustic  boundary-layer  produces  a  slight  increase 
in  the  mean  regression  rate  (1.8  percent)  on  which  the  sinusoidal  response  of  the 
regression  rate  to  the  acoustic  pressure  is  superimposed.  As  the  location  is  moved 
closer  to  the  center  of  the  chamber,  nearer  the  pressure  node,  the  acoustically 
produced  turbulence  begins  to  dominate  the  dynamics  of  the  combustion  zone.  At 
x/L=0.31 25  the  sinusoidal  behavior  of  the  response  has  nearly  vanished  as  the 
response  becomes  increasingly  rectified.  At  x/L=0.4375  the  response  is  almost 
completely  rectified  due  to  the  turbulence,  yet  the  small  acoustic  pressure 
oscillation  which  exists  produces  discernable  variations  in  the  peaks  and  troughs  of 
the  response.  At  the  pressure  node,  x/L=0.5  (not  shown),  the  response  is  completely 
rectified  and  has  a  mean  shift  of  6.6  percent,  the  maximum  in  the  chamber. 
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Figure  26.  Combined  response  of  the  condensed- phase  regression  rate  as 
a  function  of  axial  location. 


The  response  of  t he  condensed  phase  regre.-  d  -n  r  a'<  to  \ he  acoustic  boundary 
layer  is  further  investigated  iwn<:  a.  conventional  Fa>t  ''"mer  Transform  t'KFT 
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analysis  on  6  cycles  of  oscillation  containing  600  points  of  data.  The  variation  of  the 
frequency  content  of  the  response  along  the  axis  of  the  chamber  was  calculated  for 
(=1000  Hz,  rima  =  0.05,  and  two  equivalent  sand  roughness  heights,  ks=150gm  and 
ks=100|im.  Figure  27  shows  the  normalized  frequency  content  for  ks=150pm  at 
various  locations  along  the  duct  and  includes  the  responses  plotted  in  Figure  26. 
The  results  presented  at  each  axial  location  have,  in  each  case,  been  normalized  by 
the  largest  amplitude  obtained  in  the  entire  chamber  after  removal  of  the  zero 
frequency  component  corresponding  to  the  DC-shift  in  the  mean  regression  rate. 
The  primary  response  from  x/L=0  to  x/L=0.25  is  seen  to  be  that  of  the  driving 
frequency  corresponding  to  the  acoustic  pressure  oscillation.  The  lack  of  any 
harmonic  content  in  this  region  is  due  to  the  nearly  first  order  (d>=l.l)  gas-phase 
reaction.  In  the  region  from  x/L=0.3125  to  the  center  of  the  chamber  the  frequency 
content  is  characterized  by  a  steady  decrease  in  the  amplitude  of  the  fundamental 
frequency  and  a  rapid  rise  in  the  amplitude  of  the  second  harmonic  corresponding  to 
the  rectified  turbulent  response.  There  is  a  similar  appearance,  and  growth  in  the 
amplitudes  of,  the  higher  even  harmonics  along  with  the  appearance  of  some  very 
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Figure  27.  Harmonic  content  of  the  combined  response  versus  axial 
location  for  nma  =  0.05,  f  =  1000  Hz,  and  ks  =  150pm). 
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small  odd  harmonics.  The  largest  response  is  seen  to  occur  at  the  center  of  the 

chamber  due  to  the  acoustically  induced  turbulence  and  is  approximately  double  the 

amplitude  of  the  largest  pressure  response,  which  occurs  at  x/T>=0.  Examining  the  * 

frequency  content  at  the  center  of  the  chamber  reveals  that  only  a  completely 

rectified  response  with  even  harmonics  is  produced,  as  was  shown  in  the  previous 

section,  by  the  turbulent  acoustic  velocity  oscillations  alone.  It  can  be  seen  that, 

although  no  odd  harmonics  are  produced  by  either  the  pure  acoustic  pressure  or 

turbulent  acoustic  velocity  oscillations  alone,  in  combination  they  do  produce 

detectable  odd  harmonics.  The  calculated  FFT  results  corresponding  to  the  100pm 

equivalent  sand  roughness  height  were  qualitatively  identical  to  the  results 

piesemea  in  Figure  27,  differing  only  in  the  magnitude  of  the  turbulent  response. 

The  results  presented  in  Figs.  26  and  27  may  also  be  examined  to  assess 
whether  the  acoustic  oscillations  within  the  chamber  could  be  damped  or  reinforced 
by  the  propellant  response.  The  frequency  response  presented  in  Figure  27 
indicates  that  an  acoustic  oscillation  at  the  fundamental  frequency  is  capable  of 
introducing  energy  to  acoustic  modes  at  the  higher  harmonic  frequencies.  The 
response  shown  in  Figure  26  indicates  that  the  energy  added  over  one  cycle  to  the 
acoustic  oscillations  would  be  negligible  towards  the  center  of  the  chamber  due  to 
the  nearly  complete  rectification  of  the  combined  response.  The  potential  for 
combustion  instability  arising  in  an  actual  rocket  chamber,  due  to  the  combined 
response,  will  depend  upon  the  extent  to  which  the  turbulent  component  is  biased, 
and  shifted  in  phase,  with  respect  to  the  local  pressure  oscillation  bv  the  interaction 
of  the  developing  mean  axial  flow  with  the  acoustic  boundary-layer  and  turbulence 
fields.  The  relative  magnitude  of  the  turbulent  and  pressure  components  of  the 
combined  response  indicate  that  ?  h-  turbulent  component  could  potentially 
contribute  a  larger  portion  of  acou-tu  em-rgy  to  the  owraH  (lowfield  than  that  due 
to  pressure  coupling  alone 

COMPARISON  OF  ACOUSTIC  AND  STEADY-STATE  EROSIVE  BURNING 

The  experiments  of  Frump  and  l>r;ce!;’  '  <  using  center  vented  T-burners  to 
investigate  the  combustion  of  solid  propellants  in  m  oscillatory  environment, 
produced  valuable  data  concerrimg  'he  .of?  '  *  >f  :•<  m-n  erosion  throughout  the 
chambers  of  such  devices.  Kef  If.  pr.-mis  flora  <.u  >  he  .  o  *  o  o  i<>n  of  acoustic  erosion 
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throughout  the  length  of  the  chamber  of  the  T-burner.  The  results  all  exhibit  a 
maximum  increase  in  mean  regression  rate  at  the  pressure  node,  corresponding  to 
the  maximum  in  acoustic  oscillations,  and  a  minimum  in  mean  regression  rate  at 
each  end  of  the  chamber.  They  inferred  that  the  reduction  in  regression  rate  below 
its  normal  value  (often  referred  to  as  “negative  erosion”)  was  produced  by  a  negative 
pressure  response  near  each  pressure  antinode  and  a  positive  velocity  response  near 
the  pressure  node.  The  Soviet  experimenters,  Medvedev  and  Revyagin14,  verified 
the  results  of  Crump  and  Price  with  T-burner  experiments  of  their  own  and 
suggested  that  the  large  increase  in  the  mean  regression  rate  may  be  due  to  the 
turbularization  of  the  Stokes  layer  close  to  the  surface  of  the  burning  propellant. 
Figure  28  shows  the  variation  of  the  acoustic  erosion  as  a  function  of  axial  location 
within  the  chamber  for  calculations  with  ks=100|im,  125pm,  and  150pm.  All  of  the 
calculations  were  performed  at  f=1000  Hz  and  nma  =  0.05.  An  equivalent  sand 
roughness  height  of  150pm  is  found  to  produce  an  acoustic  erosion  of  approximately 
1  percent  from  very  near  the  end  of  the  chamber  to  x/L=0.3125,  where  the  erosion 
begins  to  rise  rapidly  to  a  maximum  of  6.6  percent  at  x/L=0.5.  Values  of  100pm  and 
1  25pm  equivalent  sand  roughness  heights  produce  similar  variations  of  erosion 
throughout  the  chamber,  differing  substantially  only  in  their  respective 
magnitudes.  The  acoustic  erosion  remains  zero  to  approximately  x/L=0.25  and  then 
increases  to  the  maximum  values  of  1 .6  percent  and  4  percent  at  x/L=0.5.  The  lack 
of  negative  erosion  near  the  pressure  antinodes  in  the  present  calculations  may  be  a 
result  of  the  choice  of  propellant  thermophysical  parameters  or  the  simplifying 
assumptions  made  for  the  combustion  process. 

Figure  29  shows  additional  experimental  results  concerning  erosive  burning 
obtained  by  Crump  and  Price  and  Medvedev  and  Revyagin  as  presented  in  Ref.  14. 
The  mean  shift  in  tin*  condensed  phase  regression  rate  is  plotted  as  a  function  of  (be 
velocity,  u,.,  for  both  acoustic  and  steady-state  erosive  conditions.  Here,  u0 
represents  the  centerline  velocity  of  the  steady-state  experiments  and  the  amplitude 
of  the  acoustic  velocity  oscillation  in  the  oscillatory  experiments.  A  different 
propellant  was  used  by  each  set  of  experimenters,  but  the  same  trends  were 
discerned.  First,  the  acoustic  erosion  threshold  velocity  was  found  to  be 
approximately  half  of  the  threshold  velocity  for  steady-state  erosion.  Second,  the 
acoustic  erosion  was  noted  to  be  approximately  linearly  related  to  the  maximum 
amplitude  of  the  acoustic  velocity  oscillation. 
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Figure  28. 


Figure  29. 
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Figure  30  shows  results  of  calculations  of  acoustic  and  steady-state  erosion 
for  the  propellant  considered  in  the  present  investigation,  assuming  two  values  of 
equivalent  sand  roughness  height,  ks=150pm  and  100p.m.  The  acoustic  calculations 
were  performed  at  the  pressure  node  to  insure  that  the  response  was  due  solely  to 
the  acoustically  induced  turbulence.  The  maximum  acoustic  velocity  amplitude  at 
which  a  calculation  can  be  performed  for  is  limited  by  the  corresponding  maximum 
allowable  Reynolds  number,  based  on  wall  shear  stress  and  equivalent  sand 
roughness  height,  above  which  the  validity  of  the  surface  roughness  model  is 
questionable26.  The  steady-state  calculations  were  performed  using  the  SPEC24 
code  at  the  same  operating  conditions,  and  for  the  same  propellant  formulation,  as 
the  acoustic  calculations.  The  present  calculat’  ns  predict  the  acoustic  threshold 
velocity  to  be  less  than  half  the  threshold  velocity  of  steady-state  erosion  and  the 
acoustic  response  is  seen  to  vary  nearly  linearly  with  velocity  after  the  initial 
response  appears.  Both  of  these  results  agree  well  with  the  experimental  results 
shown  in  Figure  29.  It  should  be  noted,  however,  that  the  calculated  relative  slopes 
of  the  acoustic  and  erosive  burning  curves  do  not  correspond  to  the  experimentally 
obtained  results. 
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Figure  30.  Comparison  of  calculated  acoustic  and  steady-state  erosion. 


The  ability  of  the  acoustically  induced  turbulence  to  interact  with  the 
combustion  process  is  directly  related  to  the  height  of  the  acoustic  boundary-layer, 
8aa,  which  is  inversely  related  to  the  frequency  of  the  acoustic  oscillations. 
Medvedev  and  Revyagin  reported  that  the  magnitude  of  the  acoustic  erosion  was  a 
weak  function  of  acoustic  frequency  and  suggested  it  may  be  due  to  the  decrease  in 
the  height  of  the  Stokes  layer.  Figure  31  shows  the  results  of  calculations 
performed  to  investigate  the  frequency  dependence  of  the  acoustic  erosion  for  the 
present  propellant  combustion  parameters.  The  calculations  were  performed  at 
f=1000  Hz  and  2000  Hz  with  an  equivalent  sand  roughness  height  of  ks=100pm. 
The  frequency  dependence  is  predicted  to  be  quite  strong  for  the  propellant 
thermophysical  parameters  of  tins  investigation.  For  example,  at  u(,=60  m/s,  the 
acoustic  erosion  is  predicted  to  be  approximately  9  percent  at  f=1000  Hz  and  21 
percent  at  f=2000  Hz. 
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Figure  31.  Comparison  of  calculated  acoustic  erosion  at  acoustic  oscil¬ 
lations  of  lOOt)  Hz  and  2000  Hz. 
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CONCLUSIONS 


An  analysis  of  transitional  and  turbulent  acoustic  boundary  layers  in  the 
presence  of  strong  injection  has  been  presented.  The  problem  was  approached  by 
analyzing  the  behavior  of  a  second-order  turbulence  model  rather  than  the 
traditional  Orr-Summerfeld  linearization.  An  approximate,  order-of-magnitude 
analysis  provides  a  simple  and  functionally  correct  estimate  of  transition  for  non- 
injected  acoustic  boundary-layers.  This  technique  may  prove  useful  for  estimating 
the  stability  characteristics  of  other  types  of  flows. 

The  computed  results  for  Stokes  layers  with  injection  indicate  a  substantial 
increase  in  the  acoustic  boundary  layer  thickness  for  strongly  injected  laminar  or 
turbulent  flows.  For  injected  fully  turbulent  flows,  a  pronounced  Richardson  effect 
is  obtained  and  is  accompanied  by  very  large  maximum  levels  of  turbulence  (>  20%). 
Both  axial  velocity  and  turbulence  profiles  exhibit  a  convected  wave  train  shape. 
The  turbulence  development  was  also  found  to  lag  the  rectified  acoustic  velocity  by  a 
few  degrees  at  lower  frequencies. 

Calculated  transition  results  for  conventional  Stokes  layers  are  in  agreement 
with  the  approximate  analysis  and  data  trends,  although  at  the  higher  acoustic 
amplitudes  shown  by  the  scaling  relation,  additional  nonlinear  effects  would  be 
important.  The  effect  of  injection  is  to  (nonmonotonically)  decrease  the  critical 
acoustic  Mach  number  for  transition  by  up  to  a  factor  of  about  three  for  smooth 
surfaces.  This  effect  is  frequency  dependent,  but  is  expressible  in  terms  of  the 
acoustic  and  injection  Reynolds  numbers  for  the  problem.  The  results  emphasize 
the  importance  of  the  acoustic  and  injection  Reynolds  numbers  as  the  principal 
similitude  parameters  for  near-surface  phenomena. 

A  further  appreciable  reduction  in  critical  Mach  number  (or  Reynolds 
number)  is  predicted  for  injection  velocities  with  finite,  continuous  disturbance 
levels.  This  effect  is  strong  enough  to  indicate  that  at  the  minimum  critical  acoustic 
Reynolds  number,  rather  modest  levels  of  acoustic  pressure  ratio  (~  1%)  can  induce 
significant  turbulence  levels  near  the  surface.  It  is  noted  that  the  theoretical 
predictions  concerning  acoustic  boundary-layer  transition  have  been  recently 
confirmed  by  the  experiments  of  Ma,  Van  Moorhem  and  Shortfall28,  although  at 
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relatively  low  values  of  injection  Reynolds  number.  Further  experiments  at  values 
of  Res  ~  1  are  necessary  to  evaluate  theoretical  results  in  regimes  characteristic  of 
actual  solid  rocket  combustion. 

The  potential  for  turbularization  of  a  near-surface  reaction  zone  in  ducted 
flows  such  as  rocket  chambers  therefore  theoretically  exists.  However,  whether  this 
mechanism  can  produce  instability  in  such  systems  depends  upon  several  effects  not 
considered  in  this  investigation.  For  example,  combustion  processes  and 
(acoustically)  nonlinear  behavior  such  as  thermoacoustic  streaming  could  influence 
the  nature  of  the  overall  response,  while  the  axial  mean  flow  could  bias  the 
transitional  characteristics  over  the  acoustic  velocity  and  pressure  cycles.  Further 
research  on  these  topics  is  in  progress. 

The  analysis  presented  in  this  vrork  has  attempted  to  identify  and  quantify 
some  aspects  of  the  role  of  acoustically  induced  turbulence  in  the  combustion 
response  of  solid  propellants.  Despite  several  simplifying  assumptions,  a  complex 
nonlinear  and  spatially  dependent  interaction  occurs  in  the  reaction  zone  which 
produces  a  correspondingly  complex  response  function.  Some  features  of  the 
response  have  long  been  hypothesized  (see,  for  example.  Ref  1 .)  and  are  presently 
analytically  confirmed,  viz:  “full  wave”  rectification  relative  to  the  harmonic 
acoustic  velocity  outside  the  boundary  layer,  the  occurrence  of  a  threshold  acoustic 
amplitude  for  significant  response,  and  the  appearance  of  a  time-mean 
augmentation  (“D  C  shift”)  of  the  propellant  regression  rate.  Additionally,  the 
magnitude  of  the  response  is  quite  sensitive  to  propellant  "formulation  variables”; 
surface  roughness  and  gas  phase  heat  of  reaction  (flame-zone  thickness)  being  those 
assessed  in  this  study. 

The  results  also  indicate  differences  with  some  prevailing  conceptual  aspects 
of  velocity  coupled  response.  These  include: 

1)  The  response  to  acoustic  velocity  oscillations  can  be  independent  of 
steady  state  erosive  burning,  as  demonstrated  by  the  neglect  of  an 
appreciable  axial  mean  flow  in  this  analysis.  Consequently, 
contemporary  treatments  of  velocity  response  functions  which  simply 
perturb  an  erosive  burning  response  would  omit  an  important 
additional  contribution  and  physii  al  me'  hamsm. 
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2)  For  nominal  combustion  parameters,  the  susceptibility  to  acoustic 
transition  and  propellant  response  increases  with  increasing  mean 
pressure. 

3)  The  appreciable  phase  shift  of  the  regression  rate  response  relative  to 
the  harmonic  external  axial  velocity  can  yield  a  significant  “apparent” 
pressure  response,  if  an  attempt  was  made  to  analyze  this  with 
conventional  linear  theory. 

4)  The  behavior  of  the  velocity  response,  even  at  moderate  acoustic 
pressure  amplitudes  of  a  few  percent,  is  highly  nonlinear  and  possesses 
significant  harmonic  content. 

The  present  results  have  been  obtained  utilizing  several  simplifying 
assumptions,  some  of  which  pose  significant  research  issues  in  their  own  right  (for 
examp’ e,  the  potential  for  direct  acoustic-turbulence  interaction,  or  the  potential  for 
extreme  combustion-turbulence  interaction  from  composite-propellent  diffusion 
flamelets).  One  of  the  more  intriguing  questions  posed  by  the  present  results  is  the 
influence  of  an  axial  mean  flow  interacting  with  the  acoustic  boundary-layer  and 
turbulence  fields.  In  this  case  a  substantial  biasing  of  the  propellant  response  over 
a  cycle  is  anticipated  (a  necessary  condition  for  instability).  Considering  the 
sensitivity  and  nonlinearity  of  the  acoustic  transition  process,  this  area  of  research 
could  provide  sp-  ;al  insight  into  the  coupling  mechanisms  of  producing  instability 
in  actual  devices. 

The  results  presented  in  this  work  provide  an  attempt  at  evaluating  the 
combined  response,  due  to  simultaneous  acoustic  pressure  and  (turbulent)  velocity 
oscillations,  of  the  combustion  of  a  solid  propellant.  An  order  of  magnitude  analysis 
performed  on  the  mean  enthalpy  equation  demonstrated  that  it  and  the  mean 
species  equation  remain  similar  in  the  presence  of  acoustic  oscillations  in  the  gas- 
phase  reaction  zone.  This  result  should  prove  useful  in  future  investigations  of 
similar  unsteady  problems  involving  premixed  combustion,  since  chemical  species 
transport  equations  do  not  have  to  be  solved  computationally. 
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Several  of  the  characteristics  of  acoustic  erosion  which  have  been  determined 
or  inferred  experimentally  using  T-burners  are  presently  verified.  These  include: 

1)  The  maximum  of  acoustic  erosion  occurs  at  the  velocity  antinode 
(pressure  node)  and  the  minimum  occurs  at  the  ends  of  the  chamber 
(pressure  antinode).  The  calculations  did  not  show'  the  negative 
acoustic  erosion  (pressure  response)  near  the  ends  of  the  chamber  as 
was  experimentally  reported.  The  authors  believe  this  discrepency  to 
be  the  result  of  a  simplified  modeling  of  the  propellant  combustion 
process. 

2)  Under  nominal  conditions,  the  threshold  velocity  of  acoustic  erosion  is 
less  than  half  the  threshold  of  steady-state  erosion. 

3)  Above  the  threshold  value,  acoustic  erosion  is  a  nearly  linear  function 
of  the  amplitude  of  the  acoustic  velocity  oscillation. 

4)  The  magnitude  of  the  acoustic  erosion  is  a  function  of  the  frequency  of 
the  acoustic  oscillation. 

The  results  have  also  provided  the  variation  of  several  characteristics  of  the 
combined  response  as  a  function  of  avial  location  within  the  chamber.  The 
response  of  the  propellant  was  found  to  change  from  the  purely  sinusoidal  response 
at  the  chamber  ends  to  the  completely  rectified  response  containing  only  even 
harmonics  at  the  center  of  the  chamber  Between  those  extremes,  the  frequency 
content  is  composed  of  the  fundamental  and  both  even  and  odd  harmonics,  which 
could  lead  to  the  excitation  of  other  modes  of  acon-tir  n>nilat  ion  within  the  chamber 
of  an  actual  operating  system. 
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APPENDIX  A:  TABLE  OF  ASSUMPTIONS 


Summary  of  Assumptions  Concerning  the  Gas  and  Condensed  Phases 

Gas  Phase  Assumptions 

a)  External  body  forces  and  the  coefficient  of  bulk  viscosity  are  negligible. 

b)  Effects  of  radiation  are  negligible. 

c)  Species  diffusion  due  to  thermal  and  pressure  gradient  effects  is 
negligible,  and  all  binary  diffusion  coefficients  arc  equal 

d)  The  Lewis  number  is  unity. 

e)  Specific  heats  of  the  chemical  species  are  equal  and  independent  of 
temperature. 

f)  The  mean  flow  Mach  number,  <M>,  is  negligible. 

g)  The  acoustic  boundary  layer  thickness,  5a,  is  much  smaller  than  the 
duct  radius  or  transverse  dimension. 

h)  Combustion  of  a  homogeneous  reactant  mixture  proceeds  through  a 
single  step,  irreversible,  chemical  reaction. 

i)  The  gas  phase  species-average  molecular  weight  is  equal  to  the 
molecular  weight  of  the  mixture  at  equilibrium. 

h)  The  turbulent  enthalpy  correlations  above  the  flame  zone  are 
negligible. 

Condensed  Phase  Assumptions 

a)  The  condensed  phase  is  homogeneous. 

b)  The  condensed  phase  is  stationary  tu  the  *  •  feienro  coordinate  system. 

c)  Material  proper? it  s  are  constant. 

d)  Species  diffusion  >•;  m  eligible 

e)  Radiation  efteets  at  i  h<>  mi<  •  f  ev  and  m  :  1’.*.  •ttitdr-n<e<l  phase  are 
negligible 

f)  The  surface  regression  rate  L  small  compared  to  the  gas  phase 
injection  \tT>aty. 

Subsurface  reactions  an-  not  ptesent. 

Surface  regression  mechanism  is  governed  b\  the  Arrhenius  pyrolysis 
law. 
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APPENDIX  B:  TABLE  OF  PROPELLANT  PARAMETER 


Propellant  Thermochemical  Parameters 


Tf  =  2976  K 

cp  =  1 .92  x  103  J/kg-K 

Wav  =  25.8  kg/kg-mole 

TAg  =  2.0  x  104  K 

Bg  =  6.5  x  109  kgi-*  nA  1  s  i  K  Pe 

O  =  1.10 

pg  =  0 

p  =  4.42  x  10~7  (T)0-65  kg/m/s/K0-65 
Pr  =  1 
Tn  =  300  K 

L°  =  —4.184  x  105  J/kg  (-100  cal/gm) 

Pti  =  1700  kg/m3 

c*  =  1.46xl03 

An  =  2.5  x  107  kg/m2-s 

Tax  =  1.0  x  104  K 

r’  =9.2x10-3  m/s 

T*  =  700  K 

p*  =  68  atm 

n  =  0.55 
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